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Abstract: Wc discuss the nature of the finite-temperature transition in QCD 
with Nf massless flavors. Universahty arguments show that a continuous (second- 
order) transition must be related to a three-dimensional universality class char- 
acterized by a complex Nf x Nf matrix order parameter and by the symmetry- 
breaking pattern [SU{Nf)L®SU{Nf)R]/Z{Nf)v ^ SU {N f)v mN f)v , or [C/(A^/)l® 
U{Nj)r\/U{X)v U{Nf)v/U{l)v if the U{1)a symmetry is effectively restored at 
Tc. The existence of any of these universality classes requires the presence of a sta- 
ble fixed point in the corresponding three-dimensional <l>^ theory with the expected 
symmetry-breaking pattern. Otherwise, the transition is of first order. 
In order to search for stable fixed points in these theories, we exploit a three- 
dimensional perturbative approach in which physical quantities are expanded in pow- 
ers of appropriate renormalized quartic couplings. We compute the corresponding 
Callan-Symanzik /3-functions to six loops. We also determine the large-order behav- 
ior to further constrain the analysis. No stable fixed point is found, except for Nf — 2, 
corresponding to the symmetry-breaking pattern [SU{2)l S'f/(2)^]/Z(2)y — > 
S'f/(2)v//Z(2)v equivalent to S0{4) — * 5*0(3). Our results confirm and put on a 
firmer ground earlier analyses performed close to four dimensions, based on first- 
order calculations in the framework of the e = A — d expansion. 
These results indicate that the finite-temperature phase transition in QCD is of first 
order for A^/ > 3. A continuous transition is allowed only for Nf — 2. But, since the 
theory with symmetry-breaking pattern [U{2)l <8) C/(2)jj] /[/(!) — > C/(2)y/C/(l) does 
not have stable fixed points, the transition can be continuous only if the effective 
breaking of the U{1)a symmetry is sufficiently large. 
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1. Introduction 

The thermodynamics of quarks and gluons described by QCD is characterized by 
a transition from a low-temperature hadronic phase, in which chiral symmetry is 
broken, to a high-temperature phase with deconfined quarks and gluons (quark-gluon 
plasma), in which chiral symmetry is restored. The main features of this transition 
depend crucially on the QCD parameters, such as the number Nf of flavors and 
the quark masses. In the light-quark regime the nature of the finite-temperature 
transition is essentially related to the restoring of the chiral symmetry |1|, ^, ^, 
For recent reviews see, e.g., Refs. [§|, §|, 0) 

In the vanishing quark-mass limit, the QCD Lagrangian is invariant under U (A^/)l 
and U{Nf)ji transformations involving the left- and right-handed quark flavors. Since 
U{N)l,r = U{1)l,r^[SU{N)/Z{N)]l,r and the group U{1)l®U{1)r is isomorphic 
to the group U{l)v ® U{1)a of vector and axial U{1) transformations, the classical 
symmetry group of the theory can be written as 



Uil)v ® U{1)a ® [SUiNj)/Z{Nf)]L [SU{Nj)/Z{Nf)]R 



(1.1) 



The vector subgroup U{l)v corresponds to the quark-number conservation, and it 
is not expected to play any role in the transition. The axial subgroup U{1)a is 
broken by the anomaly at the quantum level, reducing the relevant symmetry to [|I| 
[SU{Nf)i ® SU{Nf)R]/Z{Nf)v- At T = the symmetry is spontaneously broken 
to SU{Nf)v with Nj — 1 Goldstone particles (pions and kaons) with a nonzero 
quark condensate {ipip)- With increasing T, a phase transition characterized by the 
restoring of the chiral symmetry is expected at Tc] above the quark condensate 
vanishes. The symmetry-breaking pattern at the transition is expected to be 

[SU{Nf)L ® SU{Nf)R]/'L{Nf)v ^ SU{Nf)vl'L{Nf)v. (1.2) 

with order parameter given by the expectation value of the quark bilinear '^ij = 
ipL,i'4^R,j- In this picture the quark masses act as an external field Hij coupled to the 
order parameter. Note that, if the U{1)a symmetry is effectively restored at Tc, the 
relevant symmetry-breaking pattern would be 

[U{Nf)L ® U{Nf)R]/U{l)v U{Nf)y/U{l)v . (1.3) 

The effects of the anomaly tend to be suppressed in the large-Nc limit, where Nc is 
the number of colors. The symmetry-breaking pattern ( |1.3|) should be recovered in 
the limit of an infinite number of colors. 

In the case of a continuous transition, thermodynamic quantities are analytic 
functions of T as long as the external field coupled to the order parameter is nonva- 
nishing. Therefore, if the chiral-symmetry restoring transition is continuous in the 
vanishing quark- mass limit, an analytic crossover is expected for nonvanishing quark 
masses (actually for not too large values because in the heavy-quark limit the tran- 
sition is first order 0). On the other hand, if the phase transition is of first order 
in the chiral limit, it persists also for nonvanishing masses, up to a critical surface 
in the quark-mass phase diagram where the transition becomes continuous and is 
expected to be Ising-like Outside this surface, i.e. for larger quark masses, the 
phase transition disappears and we have an analytic crossover. 

The nature of the finite-temperature transition in QCD with Nf light flavors has 
been much investigated, exploiting various approaches, including numerical simula- 
tions of lattice QCD. Our understanding is still essentially based on the universality 
arguments reported by Pisarski and Wilczek in Ref. [|T|. Their main points can be 
summarized as follows: 

(i) Let us first assume that the phase transition at is continuous (second or- 
der) for vanishing quark masses. In this case the critical behavior should be 
described by an effective three-dimensional {3-d) theory. Indeed, the length 
scale of the critical modes diverges approaching T^, becoming eventually much 
larger than l/T^ that is the size of the euclidean "temporal" dimension at Tc. 
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Therefore, the asymptotic critical behavior must be associated with a 3-d uni- 
versahty class characterized by a complex Nf x Nf matrix order parameter, 
corresponding to the bilinear \E'jj = ipL,i'4'R,ji ^i-nd by the symmetry-breaking 
pattern ( p. .21) , or ( |1.3| ) if the U{1)a symmetry is effectively restored at Tc- 



(ii) According to renormalization-group (RG) theory, the existence of such uni- 
versality classes can be investigated by considering the most general Landau- 
Ginzburg- Wilson (LGW) $^ theory for a complex Nf x Nf matrix field $jj with 
the desired symmetry and symmetry-breaking pattern. The LGW Lagrangian 



for the symmetry-breaking pattern ( |1.3| ) is given by 

Cu^Nf) = Tr(9^<l>t)(9^<l>) + rTr^l^^ + ^ (Tr*^*)^ + ^Tr {^^^f . (1.4) 



The chiral anomaly reduces the symmetry to ( |1.2| ) and thus new terms must 
be added. The most relevant one is proportional to (det$^ + det$) , which is 
a polynomial of order Nf in the field $. Thus, for A'^^ > 5 such a term, and 
therefore the [/(l)A-symmetry breaking, is irrelevant at the transition. Instead, 
for Nf = 3 and Nf = 4 the new term must be added to the Lagrangian (|1.4| ), 
obtaining 

J^suiNf) = Cu[Nf) + Wo (det<l>^ + det$) . (L5) 

For Nf < 2 additional terms depending on the determinant of $ should be 
added, in order to include all terms with at most four fields compatible with 
the symmetry. Nonvanishing quark masses can be accounted for by adding an 



external- field term Hij, such as 



Tr + h.c.) . (L6) 



(iii) The critical behavior at a continuous transition is described by the stable fixed 
point (FP) of the theory, which determines the universality class. The absence 
of a stable FP indicates that the phase transition is not continuous. Therefore, a 
necessary condition of consistency with the initial hypothesis that the transition 
is continuous, cf. (i), is the existence of stable FP's in the theories described 
by the Lagrangians (|1.4| ) and (|L5|) . 

(iv) In the one-flavor case, Nf = 1, the Lagrangian CuiNf) is equivalent to the 
0(2)-symmetric real Lagrangian, which has a stable FP describing the XY 
universality class. For Nf = 1 the determinant term in the Lagrangian Csu{Nf) 
is equivalent to an external field coupled to the order parameter, which smooths 
out the singularity of the continuous transition. 

In the two-flavor case, Nf = 2, the symmetry-breaking pattern ( p..2| ) is equiva- 
lent to 5*0(4) —>■ S0{3), which is the symmetry-breaking pattern of the stan- 
dard 5'0(4)-symmetric theory for a four-dimensional vector fleld. As is well 
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known, such a theory has a stable FP, corresponding to the SO (4) univer- 
sahty class. Thus, for Nf = 2, there exists a three-dimensional universality 
class with the desired symmetry-breaking pattern. This fact tells us that, if 
the phase transition is continuous, it belongs to the SO (4) universahty class, 
i.e. it has the same critical exponents, scaling functions, critical equation of 
state, etc.... However, these arguments do not exclude that the transition is a 
first-order one. 



(v) The problem is more complex in the other cases, i.e. Nj > 3 for the symmetry 
breaking (|1.2|) and Nf > 2 for the symmetry breaking ( |1.3|) . In order to inves- 
tigate the RG flow of the Lagrangians (|1.4|) and (|1.5|) , Pisarski and Wilczek [|l| 
performed a first-order perturbative calculation within the e-expansion frame- 
work, i.e. an expansion in powers of e = 4 — d. Within this approximation, 
they did not find stable FP's. Then, extending this result to three dimensions, 
i.e. to e = 1, they argued that the finite-temperature transition in QCD is 
always a first-order one, with the only possible exception of Nf = 2. 



The e expansion provides useful indications for the behavior of the RG fiow 
in three dimensions. But the validity of the extrapolation to e = 1 of the results 
obtained near four dimensions is not guaranteed, not even at a quahtative level, 
especially if they are obtained from low-order computations. The location and the 
stability of the FP's may drastically change approaching d = 3, and new FP's, not 
present for d ^ 4, may appear in three dimensions. In several physically interest- 
ing cases low-order (and in some cases also high-order) e-expansion calculations fail 
to provide the correct physical picture. We mention the 3-d 0{N)(B0{M) invari- 
ant $^ theory describing the multicritical behavior arising from the competition of 
distinct order parameters |10, 11 the 3-d 0(A^)®0(M) $^ theory describing the 



critical behavior of frustrated spin models with noncoUinear order [T^, |T^,^ and the 
Ginzburg-Landau model of superconductors, in which a complex scalar field couples 



to a gauge field [|T^, |1^. For example, in the latter case one-loop e-expansion calcula- 
tions indicate that no stable FP exists unless the number N of real components 
of the scalar field is larger than N^. = 365. This number is much larger than the phys- 
ical value N = 2. Consequently, a first-order transition was always expected [|T^]. 
Later, exploiting three-dimensional theoretical approaches (see, e.g., Ref. JlSl) and 
Monte Carlo simulations (see, e.g., Ref. |jl6|), it was realized that three-dimensional 
systems described by the Ginzburg-Landau model can also undergo a continuous 
transition — this implies the presence of a stable FP in the 3-d Ginzburg-Landau 



^Particularly important cases are iV = 1, M = 2, relevant for anisotropic antiferromagnets, and 
M = 2, N = 3 relevant for high- Tc superconductors. 

^These models are physically realized by stacked triangular antiferromagnets, in the case N = 2, 
M = 2,3. 
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theory — in agreement with experiments |17]. Therefore, a three-dimensional analy- 
sis, not relying on extrapolations from e ^ 1 to e = 1, is called for in order to check 
the picture reported at point (v) above. 

In this paper we return to point (v), i.e. to the issue concerning the existence of 
three-dimensional universality classes with the symmetry-breaking pattern of QCD, 
and therefore the existence of stable FP's in the theories defined by the Lagrangians 
( |1.4[ ) and ( [1.51 ). In order to obtain a reliable picture of the RG flow, we perform 
a calculation in a perturbative framework in which the dimension d is fixed to the 
physical value, i.e. d = 3, and the expansion is performed in powers of appropriate 
renormalized quartic couplings. We compute the corresponding Callan-Symanzik (3- 
functions to six loops. Moreover, we study the general properties of the expansions: 
we determine the region in which they are Borel summable and compute their large- 
order behavior to further constrain the analysis. This allows us to reconstruct the Z-d 
RG flow, and in particular to check the existence of FP's, with a quite high confidence. 
We show that no stable FP exists in all cases considered at point (v) above, i.e. 
Nf > 3 for the symmetry breaking ( |1.2| ) and Nf > 2 for the symmetry breaking 
( |1.3|) . Therefore, we fully confirm the conclusions of Ref. Concerning QCD, 
these results indicate that the finite-temperature phase transition is of first order for 
Nf > 3. A continuous transition is allowed only for Nf = 2. But, since no stable FP 
is associated with the symmetry breaking [U{2)l <S) U{2)r]/U{1)v — > U{2)v/U{l)v, 
the transition can be continuous only if the effective breaking of the U{1)a symmetry 
is sufficiently large. Otherwise, it is first order. For a specific effective breaking of the 
U{1)a symmetry, we may also have a mean-field critical behavior with logarithmic 
corrections, which is associated with the tricritical point that separates the first-order 
and the second-order transition lines in the T — g plane, where g parametrizes the 
effective breaking of the U{1)a symmetry. 

Direct information on the QCD phase transitions has been provided by numerical 
Monte Carlo studies within the lattice formulation of QCD, see, e.g., the recent review 

The Monte Carlo results obtained so far are consistent with and substantially 
support the above-reported three-dimensional RG results. In the case of three-flavor 
QCD the evidence that the transition is flrst order is rather clear, see, e.g., Refs. 

0151 EDI, PII, 521. Moreover, it has been verifled that the flrst-order transition 



persists for nonvanishing quark masses up to Ising-like endpoints [0. In the two- 
flavor case recent lattice simulations have provided a rather convincing evidence of 



a continuous transition at ~ 172 MeV, see, e.g., Refs. |23|, g6|, whose 

scaling properties are substantially consistent with those predicted by the ?>-d SO (4) 
universality. However, a conclusive evidence in favor of an SO (4) scaling behavior 
in the continuum limit has not been achieved yet. Apparently, a mean-fleld critical 
behavior has not been ruled out yet. 

The paper is organized as follows. In Sec. |^ we present the general three- 
dimensional perturbative framework. In Sec. we report our calculations for a 
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theory with symmetry [U{N)l ® U{N)r]/U{1)v, i-e. for the Lagrangian ([T^).^ 
The analysis of the perturbative expansions, using also information on their large- 
order behavior, shows no evidence of stable FP's, for any A^. In Sec. ^ we consider 
the more general Lagrangian (|1.5| ), which is symmetric under the smaller group 
[SU (N) L®SU {N)r]/Z{N)v. We present a six-loop analysis for = 4, which is the 
only case for which no convincing arguments exist in favor or against the existence of 
a stable FP. Again, our analysis does not find evidence of a stable PP. Finally, in the 
Appendix we discuss in detail the phase diagram of the general $^ model realizing 
the symmetry-breaking pattern (|1.2|) for N = 2 with a complex 2x2 matrix order 
parameter, which is relevant for the finite-temperature transition in two-flavor QCD. 



2. Perturbative expansion in Landau-Ginzburg- Wilson the- 
ories 

2.1 General $^ theories 

In the RG approach to critical phenomena, many phase transitions can be investi- 
gated by considering effective LGW theories, containing up to fourth-order powers 
of the field components. The general LGW Lagrangian for an A^-component field (pi 
can be written as 

-^LGW = j ^(d^cpif + ^ + ^ J^Mijfcz (p^jMi , (2.1) 

i i ijkl 

where the number of independent parameters and Uijki depends on the symmetry 
group of the theory. An interesting class of models are those in which 0^ is the 
only quadratic invariant polynomial. In this case, all rj are equal, = r, and Uijki 
satisfies the trace condition [p^ 



^Uiiki(x5ki. (2.2) 

i 

In these models, criticality is driven by tuning the single parameter r. Therefore, 
they describe critical phenomena characterized by one (parity-symmetric) relevant 
parameter, which often corresponds to the temperature. Of course, there is also (at 
least one) parity-odd relevant parameter that corresponds to a term Yli ^i't'i that 
can be added to the Hamiltonian (|2.1| ). For symmetry reasons, criticality occurs for 
/ij — * 0. The simplest example of this class of models is the O(A^) -symmetric 
theory, which has only one quartic parameter u: 

Uijki = 2 {Sij6ki + SikSji + 6ii6jk) . (2.3) 



^For simplicity, from now on we simply write N instead of Nf. 



- 6 - 



This model describes several interesting universality classes: the Ising one for = 1 
(e.g., liquid- vapor transition, transitions in uniaxial magnetic systems, etc.), the 
XY one for N = 2 (e.g., superfluid transition in ^He, transitions in magnets with 
easy-plane anisotropy, etc.), the Heisenberg one for = 3 (e.g.. Curie transition 
in isotropic magnets); for ^ it describes the behavior of dilute homopolymers 
for large degree of polymerization. See, e.g., Refs. for recent reviews. But 



there are also other physically interesting transitions that are characterized by more 
complex symmetries: for instance, the magnetic transitions in randomly dilute mag- 
netic systems, in stacked triangular antiferromagnets, in Mott insulators, and several 



structural phase transitions. See, e.g., Refs. [^, 13^, Q for reviews. In this case, 
the corresponding LGW Lagrangian contains two or more quartic couplings, but still 
only one quadratic invariant. More general LGW Lagrangians, that allow for the 
presence of independent quadratic parameters r^, have been considered to study the 
multicritical behavior that arises from the competition of distinct types of ordering, 
see, e.g., Refs. [|Tg, 

In the field-theory (FT) approach the RG flow is determined by a set of RG 
equations for the correlation functions of the order parameter. In the case of a 
continuous transition, the critical behavior is determined by the stable FP of the 
theory, which characterizes a universality class. The absence of a stable FP is usually 
considered as an indication that the phase transition is of first order, even in those 
cases in which the mean-field approximation predicts a continuous transition. But, 
even in the presence of a stable FP, a first-order transition may still occur for systems 
that are outside its attraction domain. 

The RG flow can be studied by perturbative methods, for example by performing 
an expansion in powers of e = 4 — or a fixed-dimension (FD) expansion in pow- 
ers of appropriate zero-momentum renormalized quartic couplings. The e expansion, 
which was introduced by Wilson and Fisher is based on the observation that, for 
d = 4, the theory is essentially Gaussian. One considers the standard perturbative 
expansion, and then transforms it into an expansion in powers of e = 4 — ci. Subse- 



quently, Parisi ||3^ pointed out the possibility of using perturbation theory directly 



at the physical dimensions d = 3 and d = 2 in the massive (high-temperature) phase. 
In the following we will follow the latter approach. 

2.2 The fixed-dimension expansion 

Let us consider the generic $^ Lagrangian (|2.1| ) for an iV-component real field 



with a single parameter r, i.e. rj = r for all r, so that the quartic terms satisfy 
the trace condition (p.2|) . In the framework of the FD expansion one works directly 
in the dimension of interest, d = 3 or d = 2. In this case the theory is super- 
renormalizable, since the number of primitively divergent diagrams is finite. One 
may regularize the corresponding integrals by keeping d arbitrary and performing an 
expansion m6 = 3 — doTe = 2 — d. Poles in e appear in divergent diagrams. Such 
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divergences are related to the necessity of performing an infinite renormalization 
of the parameter r appearing in the bare Lagrangian, see, e.g., the discussion in 
Ref. p3. This problem can be avoided by replacing r with the mass m defined by 



m 



9r(2)(p2 



r(2)(o) dp^ 



[2A) 



where the function T^'^\p'^) is related to the one-particle irreducible two-point func- 
tion by 



^.;{p) = S,,T^'\p'). (2.5) 

Perturbation theory in terms of m and Uijki is finite. The critical limit is obtained 
for m — > 0. To handle it, one considers appropriate RG functions. Specifically, 
one defines the renormalized four-point couplings gtjki and the field-renormalization 
constant by 



r(4) 

ijkl 



(0) 



5,,Z^^ [m^+p^ + Oip' 

—2 

9ijkl, 



(2.6) 
(2.7) 



where ri"^...,a„ are n-point one-particle irreducible correlation functions. Eqs. ( |2.6| ) 
and ( |2.7| ) relate the mass m and the renormalized quartic couplings gijki to the 
corresponding bare Lagrangian parameters r and Uijki- In addition, one introduces 



(1-2) _ 3ijZt\ where T^^'^) is the 
one-particle irreducible two-point function with an insertion of | (pf. 



the function Zt that is defined by the relation rj^j'^'^(O) — u^^^^ 



The common zeroes of the Callan-Symanzik /^-functions 

dgijki 



Pijkl{gabcd) — 



dm 



[2.8) 



"^abcd 



provide the FP's of the theory. In the case of a continuous transition, when m — > 0, 
the couplings gijki are driven towards an infrared-stable FP. The stability properties 
of the FP's are controlled by the eigenvalues uoi of the matrix 

^ ^ijkl,abcd rj V^'"^/ 
Ogabcd 

computed at the given FP: a FP is stable if all eigenvalues Ui have positive real part. 
The critical exponents are then obtained by evaluating the RG functions 



dlnZ, 



ijkl ) 



dlnm 



Vtig 



dlnZf 



ijkl) 



dlnm 



at the stable FP ^* r] = r]^{g*jki), and u = [2 - r]^{g*- 



(2.10) 
)]~^. From the 



Jijki) + ^t[gijki 

perturbative expansion of the correlation functions T^'^\ T^^\ and r(^'2) and using 
the above-reported relations, one can derive the expansion of the RG functions Pijki, 
7]^, and r]t in powers of gijki- 
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2.3 Resummation of the series 

FT perturbative expansions are divergent. Thus, in order to obtain accurate results, 
an appropriate resummation is required. In the case of the 0(A^)-symmetric models, 
accurate results have been obtained from the analysis of the available FD six-loop 
^ by exploiting Borel summability (proved in dimension d < 4 [^) and the 



series 



knowledge of the large-order behavior of the expansion, see, e.g., Refs. |5S[. If we 
consider a quantity S{g) that has a perturbative expansion^ 

S{g)^J2"^9\ (2.11) 

the large-order behavior of the coefficients is generally given by 

Sk ~ k\ i-a)'' e [I + 0{k-^)\ . (2.12) 

Borel summability of the series for (? > requires a > 0. The value of the constant a 
is independent of the particular quantity considered, unlike the constant h. They can 
be determined by means of a steepest-descent calculation in which the saddle point 
is a finite-energy solution (instanton) of the classical field equations with negative 



coupling see also Refs. p9|, 0]. In three dimensions and for any A^, 



0.00881962..., see, e.g., Ref. [2^. In order to resum the perturbative series, one may 



introduce the Borel-Leroy transform B{t) of S{g), 

POO 

Sig) = / t^e-*5(t), (2.13) 
Jo 

where c is an arbitrary number. Its expansion is given by 

Be.M) = V -t\ (2.14) 

The constant a that characterizes the large-order behavior of the original series is 
related to the singularity ts of the Borel transform B{t) that is nearest to the origin: 
tg = —1/a. The series i?exp(^) is convergent in the disk |t| < = 1/a of the complex 
plane, and also on the boundary if c > 6. In this domain, one can compute B{t) 
using Be^p{t). However, in order to compute the integral ( |2.13| ), one needs B(t) for 
all positive values of t. It is thus necessary to perform an analytic continuation of 
Bcxp{t). For this purpose one may use Fade approximants to the series ( p.l4|) [p5|] : 
This is the so-called Pade-Borel resummation method. A more refined procedure 
exploits the knowledge of the large-order behavior of the expansion, and in particular 
of the constant a in Eq. ( p.l2| ). One performs an Euler transformation 

, , + at — 1 , ^ 

yit) = >— (2.15) 

y'l + at + 1 



Here the coupling g is normalized according to Eq. (2.7) with gijki 
^ (SijSki + SikSji + SiiSjk) g. 

^ More precisely, the saddle-point fs{x) is a solution of the equation Afs = fs — fg- One 
show that a is related to its classical action, indeed a = — l/S*/ where Sf = — § / cfix pOt - 
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that allows to rewrite B{t) in the form 



Bit) = J2hm]'. (2.16) 

k 

If all the singularities of B{t) belong to the real interval [—oo,—tc\, the expansion 
( |2.16D converges everywhere in the complex t-plane except on the negative axis for 



t < —tc- After these transformations, one obtains a new expansion for the function 

S{g)^ / dte-H^Y.f^'^y^^st- (2.17) 
•^0 k 

This sequence of operations has transformed the original divergent series into an 
integral of a convergent one, which can then be studied numerically. However, the 
convergence of the integral (|2.17|), that is controlled by the analytic properties of 



S{g), is not guaranteed. Indeed, since generic RG quantities have a cut for g > g* 



[p^^ , |4^ , |43| , where g* is the FP, the integral does not converge for g > g*. 

The method can be extended to the case of more general LGW theories, with 
several quartic couplings, see, e.g., Refs. |^ In the case of k independent quartic 
Lagrangian terms, a generic quantity S depends on k renormalized couplings. Its 
expansion can be written as 

Sigi, ...,gk) = J2 <'n...^^9^-9i: = Y.'-9\ (2.18) 

ji,---,jk n 

where 



M'-9i\ h = 9,l9- (2.19) 

ji+---+ife=" 

The large-order behavior of the coefficients s^, and therefore the singularities of the 
Borel transform, depends on the ratios ^j, but it can be still determined by steepest- 
descent calculations. Borel summability requires that the Borel trasform does not 
have singularities on the positive real axis. Under the assumption that the group 
and the space structure of the relevant saddle points are decoupled,^ it is possible to 
show that the region where Borel summability holds is directly related to the region 
of stability of the Lagrangian potential with respect to the bare quartic couplings: 
with the normalizations implicit in Eq. ( |2.7| ), it is enough to replace the bare quartic 
couplings with their renormalized counterparts. For instance, in 0(A^)-symmetric 
models the potential is stable only for m > and therefore the expansion is Borel 
summable only for > 0. Once the large-order behavior of the coefficients Sn has 
been determined, i.e. 

s„~n![-A(^l,...,^fe)^n^ (2.20) 



^This fact has been proved to hold in 0(iV)-symmetric models 
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one may use a conformal- mapping method based on Eq. ( 2.15|) to resum the series 



in as in the 0(A^)-symmetric case. Again, A{gi, gk) does not depend on the 
quantity at hand. Alternatively, one may use the Pade-Borel method, employing 
Pade approximants to analytically extend the Borel transform. 

3. The U{N)l ® U{N)r effective theory 
3.1 The effective LGW theory 

The most general $^ Lagrangian that is symmetric under U{N)l ® U{N)fi transfor- 
mations is given by 

= Tr(9^$t)(9^$) + rTr $t$ + ^ (Tr + ^Tr ($^$)' , (3.1) 

where the fundamental field is a complex N x N matrix. There is only one 
quadratic invariant and therefore the trace condition ( |2.2| ) is satisfied by the quartic 
terms. The Lagrangian ( p.l|) is invariant under the transformations $ —>■ A^B, A, 
B G U{N), and the corresponding symmetry group is [U{N)l®U{N)ji]/U{1)v- The 
condition f o > is necessary to realize the symmetry breaking 

[U{N)l ® U{N)r]/U{1)v ^ U{N)v/U{l)v, (3.2) 

where U{N)y is the group of transformations $ — A G U{N). Indeed, 
for r < and f o > the minimum of the quartic potential, which has the form 
^ij = is symmetric under U{N) transformations, thus leading to the desired 

symmetry-breaking pattern. For vq < the minimum of the potential is obtained 
for $jj = (poSiaSja, where a is an arbitrary integer, 1 < a < A^. As one can easily 
check, this leads to a different symmetry-breaking pattern. Note that for A^ = 1 the 
two quartic terms are equal, and the Lagrangian £r/(Ar) becomes equivalent to the 
standard f/(l)-symmetric $^ Lagrangian for a complex field with coupling uq + Vq, 
which in turn is directly related to the 0(2)-symmetric real Lagrangian. 

As explained in Sec. p.2| , the theory is renormalized by introducing a set of 
zero-momentum conditions for the one-particle irreducible two-point and four-point 
correlation functions: 

^al^iP) = = <^^^<^^^^Z^' K + 0{p')] , (3.3) 



(0) = ^ 

aia2,bib2,ciC2,dl<i2^ ' r^f r-=ht 



(3.4) 



zero mom. 



b^\^a2^%^^J>^c^C2^^d 
= Z^ Tfl {uUaia2,bib2,ciC2,d-id2 ~^ 'V^aia2,bib2,ciC2,d-id2^ 

where F is the generator of the one-particle irreducible correlation functions (effective 
action), and 

Uaia2,bib2,ciC2,did2 ~ {^aiCi^bidi^a2C2^b2d2 ~^ ^aidi ^6iCi '^02^2 '^f'2C2 ) 5 
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Vaia2,bib2,ciC2,did2 2 i^aici^bidi^a2d2^b2C2 ~^ ^aidi^bici^a2C2^b2d2) ■ 

The FP's of the theory are given by the common zeroes of the Callan-Symanzik 
/3-functions 

du 

dm 



, Pv[u,v) = m— 

UO,VO 



(3.5) 

UO,VO 



Their stabihty is controlled by the eigenvalues of the matrix 



^ _ f OPu/du dPu/dv\ 

A FP is stable if all the eigenvalues of its stability matrix have positive real part. 
According to RG theory, a stable FP in the region f > would imply the existence of 
a universality class describing continuous transitions characterized by the symmetry 
breaking ( p.2| ). Therefore, in order to establish the existence of a 3-d universality 
class with symmetry-breaking pattern (|3.2|) and a complex N x N matrix order 
parameter, one must search for stable zeroes of the /3-functions. 

One can easily identify two FP's in the theory described by the Lagrangian CuiN), 
without performing any calculation. The first one is the Gaussian FP for u = v = 0, 
which is always unstable. Since for vq = the Lagrangian CuiN) becomes equivalent 
to the one of the 0(2A^^)-symmetric model, the corresponding 0(2A^^) FP must exist 
in the v = axis for m > 0. One- loop e-expansion computations do not find other 
FP's [0. The location of the 0(2A^^) FP can be obtained by using known results for 
the 0(M)-symmetric theories. Estimates relevant for our study, i.e. for 0(8), 0(18), 
0(32) and 0(72), corresponding to = 2,3,4,6, can be found in Refs. [^, ^]. 



Moreover, using the results of Refs. |TT|, ^ on the stability of the 0(iV)-symmetric 
FP under generic perturbations, one can already conclude that the 0(2iV^) FP is 
unstable for any > 2. Indeed, the term Tr ($"^$) in the Lag rangian >Cf7(jv) is a 
particular combination of quartic operators transforming as the spin-0 and spin-4 
representations of the 0(2A^^) group, and any spin-4 quartic perturbation is relevant 
at the 0(M) FP for M > 3 ]n], 0. The RG dimension 1/4,4 of the spin-4 quartic 
operators at the 0(M) FP can be computed by using the six-loop FD results and the 
five-loop e-expansion results for the related crossover exponent 04,4 = 1/4,4// reported 
in Ref. fill. We obtain the estimates ?/4,4 = 0.386(7), 0.673(6), 0.808(5), 0.914(3) 



respectively for 0(8), 0(16), 0(32), and 0(72), that are the cases interesting for this 
work. Since 1/4,4 > 0, the 0(2A^^) FP is unstable in all cases. 

3.2 Fixed-dimension expansion in three dimensions 

In order to investigate the existence of new FP's, we computed the FD expansion 
of the 3-d Callan-Symanzik /5-functions ( |3.5| ) to six loops. This is rather cumber- 
some, since approximately one thousand Feynman diagrams must be evaluated. We 



- 12 - 



employed a symbolic manipulation program, which generated the diagrams and com- 
puted the symmetry and group factors of each of them. We used the numerical results 
compiled in Ref. |^ for the integrals associated with each diagram. 
We report the series in terms of the rescaled couplings u and v, 

16vr _ 167r _ 

The six-loop expansions of the /3-functions corresponding to these rescaled variables 
are given by 

^ /- -X - -2 4Ar 3 o 190 + 82A^2 

P,{u, v) = -u + u+ ^^uv + ^^^v - (3.8) 



AOON 2 370 + 46Ar2 ^ ^ ^ 1 • ■ 



27(4 + Ar2)2 27(4 + Ar2)2 {A + mf (4 + A^2)^+i 



= + ITiV^^' + ITiV^' - 27(4 + iV2)2 ^ ^ (3.9) 



400iV __2 136 + 28A^2 „, 1 



27(4 + Ar2)2 27(4 + Ar2)2 '^{A + N^f+i 



The coefficients b^j and 6^^, with 4 < z + j < 7, are reported in the Tables |I] and ||, 
respectively. Actually, we also computed the RG functions 77,/, and rjt, cf. Eq. (|2.1CI| ), 
to six loops. But, since we will not use them, we do not report their series. Of course, 
they are available on request. 

We have done a number of checks by using known results available in the litera- 
ture. For V = 0, the expansion of jSu in powers of u reproduces that of the /3-function 
of the 0(2A^2)_gy]2i]2ietric model that can be found in Refs. 53| to six loops. 



Since we computed the series for generic U{N)l ® U{M)ji models, we also verified 
nontrivial relations holding for M = 1 and any A^, such as 

P^iz + y,z-y;N,M = l)+(3,{z + y,z-y;N,M =1) = P{z; 2iV), (3.10) 

where (3{z; 2N) is the /3-function of the 0(2A^) model and we are considering un- 
rescaled couplings normalized as in Eq. (p.7|). 

In order to resum the perturbative expansions, we use the conformal-mapping 
method described in Sec. pl3| . For this purpose we compute the function A defined 
in Eq. (|2.20|) by studying the saddle-point solutions (instantons) of the classical field 
equations. Under the assumption that the group and space structure of the saddle- 
point solutions decouple, the computation can be easily reduced to the case of the 



one-component $ theory described in detail in, e.g., Refs. ||29|, Here we only 



report the results, without giving the details of the calculations which are rather 
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standard. Consider a generic quantity S{u, v) = ^-^^ CijU^v^ and the corresponding 
expansion 

S{,gu,gv) = ^Sn{u,v)g''. (3.11) 

n 

The large-order behavior is given by 

Snr^n\[-A{u,v)fn\ (3.12) 

r ^1 

A(u, w) = AivMax u + v,u + — , (3.13) 
247r 

An = ^^a, (3.14) 
where a ^ 0.00881962. The expansion is Borel summable in the region 

M + ^>0, u + v>0, (3.15) 



which, as already mentioned in Sec. |273| , is related to the stability region of the quartic 
potential in the Lagrangian CuiN), which is given hj Uq + Vq/N > and Uq + Vq > 

i- 

3.3 Analysis of the series and results 

The knowledge of the large-order behavior of the series allows us to employ the 
conformal-mapping resummation method. As already mentioned in Sec. p.3| , this is 
essentially an extension of the method already applied in Refs. pTl BH| to resum the 



perturbative expansions in 0(A) theories. Explicitly, given a series 

i?(u,t;) = (3.16) 

hk 

we rewrite it as 

R{u,v) = R{gu,gv)\g^-^ = R{u,v;g = 1), (3.17) 

where 

R{u,v;g) = ^rr,{u,v)g'', r„(M,w)= ^Rhku''v\ (3.18) 

n h+k=n 

Then, we consider approximants E{R)p{u, v; b, a; g) defined by 

E{RUu,v;b,a;g) = y2Bkiu,v;b,a) / dtt'e~\ , (3.19) 

where 

v( \ \/^+^A{y, z)-l 

Yix; y,z) = V , (3.20) 



and A{y, z) is given in Eq. ( p.l3| ). The coefficients Bk{u, v; b, a) are determined by the 



requirement that the expansion of E{R)p{u, v; b, a; g) in powers of g gives R{u^ v; g) 
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4.0 99.8202086 + 79.8954292 + 16.4329798 + 0.674471379 

3.1 329.227699 N + 108.894438 + 6.64687835 

2.2 366.245988 + 251.386542 + 41.2000397 + 0.310944602 

1.3 280.991858 N + 83.3331935 + 3.27130728 
0,4 65.0697173 + 40.0695380 + 5.95052716 N"^ 



5,0 


-458.051683 - 415.773535 - 


1 1 1. 135357 iV*- 


8.64375912 N'^ + 0.0778229492 


4,1 


-2064.17070 N - 880.190454 


- 85.6254168 


+ 1.35288201 


3,2 


-2647.21638 - 2601.32272 - 


492.107570 TV^ - 


1.80697817 TV^ 


2,3 


-3634.53973 N - 1183.42261 A^^ 


- 70.4182766 


- 0.430339129 


1,4 


-1349.56809 - 1084.80505 - 


197.491458 A^4- 


2.65955034 A^^ 


0,5 


-336.753422 7V- 103.463574 


-4.81880464 





6,0 2596.27212 + 2603.95922 + 822.999966 + 86.8071683 + 0.911941440 



+0.0256180937/^1" 

5.1 14967.9450 N + 7570.43530 + 1000.72434 + 13.4051039 + 0.625521452 

4.2 21496.0853 + 26926.0889 + 6177.06793 + 216.477923 + 4.80379167 

3.3 44213.7572 N + 17686.3589 + 1767.24483 + 27.2537306 N'' 

2.4 20802.2607 + 22953.0190 + 5124.40912 + 171.846652 iV^ + 0.0681850017 

1.5 12267.0624 N + 4868.45580 + 459.187339 + 2.19119772 N'^ 

0,6 1359.90833 + 1466.49761 + 303.552043 A^^ + 5.48047785 

7^0 -16956.2524- 18511.7426 A^^ _ 5735 47522 _ 924.744197 iv^ - 32.1576679 A^® 
+0.314797775 N^^ + 0.0117120916 A^^^ 

6.1 -119503.736 A^ - 69794.4619 A^^ _ 11896.6226 A^^ - 458.710404 + 6.58801236 
+0.363425149 A^" 

5.2 -188784.360 - 284537.687 A^^ _ 78459.3804 A^^ - 4739.14953 A^^ + 22.5892998 A^^ 
+3.03196500 TVio 

4.3 -528347.831 A^- 257216.321 A^^ - 32795.4801 A^^ - 331.733950 A^^^ + 11.6377106 A^^ 

3.4 -292480.092 - 416299.806 A^^ „ 103998.225 A^"* - 4515.48789 A^^ + 8.83297998 A^^ 

2.5 -298671.340 A^ - 136395.073 A^^ _ 16227.9484 - 201.034856 A^^^ - 0.500032757 

1.6 -59511.6725 - 77035.0275 A^^ _ 18555.3684 A^^ _ 766.457739 A^^ - 3.10874638 A^^ 
0,7 -13984.0313 A^ - 6024.80289 A^^ _ 653.107647 A^^ - 5.22722079 A^^ 

Table 1: The coefficients bfj, cf. Eq. 



to order p. For each value of b, a, and p, E{R)p{u, v; b,a] g = 1) provides an estimate 
of R{u,v). The parameters b and a are arbitrary and can be used to optimize the 



resummation, see, e.g., Refs. ^ For example, one may require that the 
results are least dependent on the number p of terms. Usually one considers values 
in the region < 6 < 15 and — 1 < « ^ 5. The dependence of the results with 
respect to variations of b and a around their optimal values provides indications of 
the uncertainty. 

We have applied the conformal-mapping resummation method to the six-loop 
series of the /5-functions, and then determined their zeroes. The zeroes of the /3u and 
Py are shown in Fig. |l| for the cases = 2,3,4,6. In particular, P^{u,v) turns out 
to vanish only along the v = axis for any A^. The figures show that no other FP 
exists beside the Gaussian and the 0(2A^^) FP's located along the v = axis, and 
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4.0 

3.1 234.666985 + 100.520648 + 7.96126092 - 0.625553655 

2.2 426.449736 N + 97.5987129 - 2.25343030 

1.3 254.535978+ 144.854814 iV2 + 20.3052050 

0,4 66.3459284 N + 18.8271989 + 0.560179202 A^^ 



5,0 









4,1 


-1258.17306 - 606.687101 - 


72.7688516 iV*- 


1.08252821 - 0.287326259 


3,2 


-3271.88287 N - 894.130891 


- 25.2461704 


- 1.55153180 iV^ 


2,3 


-2559.77598 - 2108.80199 - 


375.743846 iV* - 


2. 13233677 A^*^ 


1,4 


-1708.58378 N - 585.048888 


- 39.4757354 TV^ 




0,5 


-266.102670 - 237.227857 7V2- 


44.6632567 iV*- 


0.496927343 TV^ 


6,0 










5.1 8093.66025 + 4354.56520 + 652.419979 + 13.7782897 - 1.54366403 
-0.159052164 iVi° 

4.2 27834.1272 N + 8658.04010 + 371.284973 - 18.1481024 N"^ - 1.18751905 

3.3 26250.5520 + 27345.7269 + 4995.12745 - 63.2324448 - 3.26781317 iV^ 

2.4 30725.1752 N + 11079.3449 + 816.715993 - 8.19919537 N'^ 

1.5 8470.20827 + 8744.03485 + 1815.09720 + 39.6191261 A^*^ 
0,6 2108.57715 A^ + 797.860452 A^^ + 68.8157914 A^^ + 0.284187609 A^^ 
7^0 

6.1 -58941.8989 - 34933.0997 A^^ _ 6302.78433 N'^ - 295.542877 A^^ + 0.957520387 A^^ 
-1.27213433 - 0.0997279331 A^^^ 

5.2 -254016.821 A^ - 89637.8297 A^^ _ 6227.84400 A^^ + 32.2770587 A^^ - 14.8328121 A^^ 
-0.951111299 A^" 

4.3 -275746.727 - 345193.386 A^^ _ 70372.4321 N'^ - 538.101215 A^^ - 67.5278644 A^*^ 
-3.69213913 Afi" 

3.4 -476545.568 A^ - 197140.916 A^^ _ 19569.5278 A^^ - 52.0633373 A^'^ - 8.74103639 A^^ 

2.5 -174726.233- 224936.606Ar2- 53193.01927V4- 2001. 32523Ar6- 7.87773734Af* 

1.6 -96977.9325 A^ - 44001.9263 - 5202.93250 A^^ _ 65.8929276 A^'^ 

0,7 -8282.26037 - 11331.8800 A^^ - 2764.77952 A^^ _ 113.637842 A^^ _ 0.318785501 A^^ 

Table 2: The coefficients 6^-, cf. Eq. 



already predicted by general considerations, see the discussion at the end of Sec. pT| . 
The M-coordinate of the 0{2N^) FP is ^^(jv^) ^ 1.30, 1.19, 1.12, 1.06 respectively for 
= 2,3,4,6 (these estimates are in agreement with those that can be found in 
Refs. |]45| , ^ for the position of the FP in 0(M)-symmetric theories). The curves 
appearing in Fig. |^ are quite robust with respect to the order of the series and to 
variations of the resummation parameters. None of the approximants of the series 
for f3y, from four to six loops, shows zeroes for \v\ > in the region in which the 
resummation is still reliable, i.e. for —2 ^ u,v 4, which is quite a large region if 
compared with the location of the 0(2A^^) FP along the v = axis. We recall that 
the relevant region is the one with vo > that corresponds to -u > 0. Concerning the 
zeroes of v), all zeroes obtained by changing the order (from four to six loops) 
and by varying b and a in the above-mentioned range, lie within the width of the full 
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lines. Moreover, perfectly consistent results have also been obtained by employing 
the Pade-Borel method, which assumes only Borel summability, without using the 
information on the position of the Borel-transform singularities. The existence of 
FP's outside the region where the resummation is reliable is quite unlikely. 

In conclusion, we do not find evidence of new FP's. According to RG the- 
ory, this means that a consistent model for a three-dimensional continuous transi- 
tion characterized by the symmetry-breaking pattern [[/{Nf)^ ® U{Nf)fi]/U{l)v — ^ 
U{Nf)v/U{l)v does not exist for any N > 2. Therefore, the phase transition in such 
systems must be a first-order one. Our results confirm and put on firmer ground ear- 
lier claims based on a first-order perturbative calculation within the e-expansion 
scheme. 



4. The SU{N)l SU{N)r effective theory 
4.1 General considerations 

If the U{1)a symmetry is explicitly broken by the anomaly, the axial U{1)a symmetry 
is reduced to Z(iV)^ If we rewrite the symmetry-breaking pattern ( |3.2| ) in the 
form U{1)a ® [SU{N)/Z{N)]l ® [SU{N)/Z{N)]r SU{N)v/Z{N)v and replace 
U{1)a with Z{N)a, we obtain 

[SU{N)l ® SU{N)r]/Z{N)v SU{N)v/Z{N)v. (4.1) 
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T T 

Figure 2: Sketches of possible pliase diagrams witli bicritical and tetracritical points, case 
(a) and (b) respectively. Full lines correspond to continuous transitions with symmetry- 
breaking pattern SO{4) SO{3), which are therefore the relevant ones for QCD. Dashed 
lines correspond to Ising transitions, not realized in QCD. Thick lines indicate first-order 
transitions. 

The LGW theory with a complex N x N matrix order parameter and such a breaking 
pattern is obtained by adding new terms to the Lagrangian Cu{n)- The most relevant 
one is proportional to the determinant of $ and thus we obtain the effective three- 
dimensional theory 

C-u{N) + Wq (det<l>t + det$) . (4.2) 

Here wq is related to the breaking of the U{1)a symmetry. 

For = 1 the added term is equivalent to an external field coupled to the order 
parameter. Therefore, it smooths out the continuous XY transition realized in its 
absence, leaving us with an analytic crossover. 

For N = 2 the phase diagram is more complex and is discussed in detail in the 
Appendix. In this case, the effective Lagrangian containing all possible interactions 
with at most four fields is given by 

Csui2) = TT{d,<l>^){d,<l>) + rTr$t$ 

+Wo (det<l>t + det$) + — (Tr$t$)2 + !^Tr 

+ ^ (Tr<l>t$) (det<l>t + det'l') + ^ [{det^^)^ + {det^f] , (4.3) 

where is a complex 2x2 matrix. 

Since the Lagrangian ( |4.3| ) has two mass terms, one expects several transition 
lines with a multicritical point. In the case relevant for QCD the multicritical point 
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should be identified with the U{2)®U (2) theory, and therefore it should correspond to 
wq = 0, xo = 0, = 0. As we discussed in the preceding Section, such a multicritical 
point corresponds to a first-order transition. Since a first-order transition is generally 
robust with respect to perturbations, it should persist even for a small breaking of the 
U{1)a symmetry. Therefore, on the basis of the analysis presented in the Appendix 
we expect two possible phase diagrams as a function of the temperature T and of 
the U{1)a breaking g, see Fig. ^. In the bicritical case, see Fig. |^ (a), the continuous 
phase transitions are associated with the symmetry-breaking pattern ( |4.1D which is 
equivalent to 5*0 (4) — > 5*0(3) 0]. Therefore, according to universality arguments, if 
two-flavor QCD undergoes a continuous transition, its critical behavior at should 
be that of the 3-d SO (4) universality class [|l], ||, ||, which has been accurately studied 
in the literature, see, e.g., the recent Refs. ^ |5^ and references therein. 

In the tetracritical case, see Fig. ^ (b), only two continuous transition lines are 
compatible with the symmetry-breaking pattern expected in QCD, those indicated 
with the full lines, see the discussion in the appendix. Note that the second-order lines 
end at particular values gtr which are tricritical points characterized by mean-field 
behavior with logarithmic corrections. In conclusion, the transition may be of first- 
order if the U{1)a breaking is small, second-order in the opposite case, and may be of 
mean-field type for a very specific value of g. Of course, a specific study of QCD is 
needed to identify which case is effectively realized. Lattice simulations of two-flavor 
QCD favor a continuous transition, see, e.g., Refs. |^ ^ |25|, |2^, with scaling 
properties that are substantially consistent with the 3-d SO (4) universality class. In 
the large- A'^c limit, A^ being the number of colors, the anomaly effects are suppressed 
by powers of l/Nc- Therefore, the finite-temperature transition is expected to be 
of first order. This picture is also supported by the fact that the finite-temperature 
transition in pure SU(A'c) gauge theories is of first order for A'c > 3 and in particular 



in the large- A^c limit, see, e.g., Refs. [^, |5^, and the quark contributions are expected 



to be suppressed by a factor 1/Ac for large Nc- This considerations suggest that for 
A"c > A^min — Monte Carlo simulations indicate A^min > 3 — the finite-temperature 
transition with two flavors changes its nature from continuous to first order. 

For A^ = 3, the determinant is cubic in the field $. Usually, a cubic term drives to 
a first-order transition, and therefore the addition of the determinant is not expected 
to soften the original first-order transition to a continuous one. Lattice simulations 



of QCD confirm this picture, see, e.g., Refs. |Tg, |T^, |20|, ^ ^2 . 

For A^ = 4 the determinant is a quartic-order term, giving rise to a generalized 
LGW theory with three quartic parameters. The effects of the determinant in 
this case will be discussed more carefully in Sec. \l.2\ As we shall see, the determinant 



^The effective J7(l)A-syniinetry breaking at finite temperature has been investigated on tlie 
lattice. The U{1)a symmetry appears not to be restored at Tc. However, the effective breaking of 
the axial U(1)a symmetry appears substantially reduced especially above Tc, as inferred from the 
difference between the correlators in the pion and 6 channels. See, e.g., Refs. [bSl M, p5l p6[. 
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does not give rise to new FP's, and therefore the phase transition remains first order 
also in the presence of a substantial effective breaking of the U{1)a symmetry. 

Finally, for N > 5 the determinant term is irrelevant because it gives rise to 
polynomial terms of degree higher than four. Therefore, for N > 5 the nature of the 



transition should not change. The lattice results for = 6 of Ref. |T9[ are consistent 
with this picture. 

4.2 Fixed dimension expansion in the case = 4 

As argued above, the theory for A^ = 4 is the only case in which it is necessary 
to investigate the effect of the determinant, since no arguments exist against the 
possibility that such a term softens the first-order transition expected in its absence. 
The three-dimensional RG flow of the Lagrangian Csu{n) for A^ = 4, cf. Eq. ( [4.2|) , can 
be studied by a straightforward extension of the calculations done for the Lagrangian 

^U{N)- 

Beside the renormalization conditions (13.31) and (B.4|), we add a further relation 



.(4) 



' (0) 

aia2,bib2,ciC2,did2\ J 



zero mom. 



(4.4) 



where e^ki is the completely antisymmetric tensor (61234 = 1). 

We computed the perturbative expansion of correlation functions T^'^\ T^^^ and 
f "^^^ to six loops. Beside the rescaled variables u and v, cf. Eq. ( |3.7|) with A^ = 4, we 
introduce w defined by = Anw/b. The corresponding /5-functions are given by 

puiu, V, w) = -u + u^ + ^uv + + - ^^u^ - ^u'^v (4.5) 
553 9 41 9 1 -■^ 2 9 V— \ -i- 



5400 225 25 75 

i+j+k>4 



P^{u, V, w) = -v + ^uv + '^v^ - - -^^u^v - ^uv^ (4-6) 
25 1350 45 ^ 

i+j+k>A 



(3n,{u,v,w) = -w + ^uw - ^vw - ^^u^w (4.7) 
11 61 o 4 



1350 2700 225 

i+j+k>4: 

The coefficients 6^;., b'"^^^, and 6^-^, with 4 < i+j + k < 7, are reported in the Table ^ 
Note that [3^ and [3^ are even functions of w, while [3u, is an odd function of w. 
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We have also determined the large-order behavior of the series. Writing 

Sign, gv, gw) = ^ v, w)g'^, (4. 



we found 



Sn-- n\[-A{u,v,w)fn\ (4.9) 

. , _ 6™ 

Aiu, V, wj = Max 

5 



_ _ _ V \w\ 
u + v,u + - + — 



The series turn out to be Borel summable for 

u + v>0, n + ^-^>0, (4.10) 

which is again the stability region of the quartic potential (with respect to the bare 
couplings Mo, fo, and wq). 

In order to resum the series of the /^-functions and search for their zeroes, we 
applied the same method employed in Sec. \i.'S\ . The analysis of the series clearly shows 
that the new /9-function jSyn vanishes only for w = 0. As a consequence, the FP's can 
only lie in the plane w = 0, and therefore they are those of the f/(4) l ® t/(4)ij theory, 
see Sec. ^A\ , which are unstable. Since no stable FP is found, the phase transition is 
expected to be first order for any value of Wq, and therefore also in the presence of a 
substantial breaking of the U{1)a_ symmetry. 

Acknowledgments 

We thank Pasquale Calabrese and Luigi Del Debbio for useful and interesting dis- 
cussions. 

A. The phase diagram of the SU{2)l SU{2)r theory 

In this appendix we discuss in more detail the phase diagram of the theory with 
Lagrangian ([4.3|). It is easy to map Csu{2) into a new one, with two 4-dimensional 
vector fields and SO (4) symmetry. If 

^ = a + il3 + i{A + iB) ■ a, (A.l) 

where a, (3, A, and B are real and a are the Pauli matrices, we define two fields 

0i = 2(a,l), 02 = 2(/3,5). (A.2) 

Then Eq. ( [4.3| ) can be rewritten as 

+ jMif + %m + ■ 02)^ (A.3) 
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(a) , (b) 




Figure 3: Sketches of the possible phase diagrams in the mean-field approximation: case 
(a) is characterized by a bicritical point, case (b) by tetracritical point. The dashed and 
thick lines represent second- and first-order transitions, respectively. 



with 

ri = r + Wo, 

r2 = r - Wo, 
3, 

ao = -(2mo + vo + 2xo + yo), 
3 

bo = -(2mo + vo- 2xo + yo), 
3 

Co = 2(2^0 + ^vo - yo), 

do = -3{vo + yo)- (A.4) 

The Lagrangian (|A.3|) is stable for > 0, 60 > 0, Cq > —2\/aobo, and Cq + do > 
—2\/aobo. It is easy to identify the symmetries of the Lagrangian (|A.3|) . It is invariant 



under 0(4) rotations of both fields 0i and 02 and under the Z(2) transformations in 
which 01 —01 and 02 "~02 independently (if the two transformations are applied 
simultaneously, we obtain an 0(4) transformation). Therefore, the symmetry group 
is Z(2) ® 0(4). It is interesting to understand these symmetries in the context of the 
original model ( |4.3| ). First, there is the invariance under $ — > A^B, A, B ^ SU{2) 
that corresponds to the symmetry [SU{2) S'f/(2)]/Z(2), which is equivalent to 
SO (4). The other two Z(2) symmetries correspond to the transformations $ — > <l>^ 
and $ — » Note that these two transformations are symmetries of the effective 
theory for any A^. We have not mentioned them before because they do not play any 
role at the QCD transition. 

The Lagrangian ( A. 3 ) has two independent quadratic (mass) terms and thus 



it describes critical and multicritical transitions. It is easy to determine the phase 
diagram of the theory in the mean-field approximation, extending the analysis of Liu 
and Fisher p3 . In this case, one assumes that the basic fields are space-independent. 
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thereby neglecting the effect of fluctuations. In this approximations two different 
phase diagrams are found, see Fig. ^ In one case, there are three phases: 1) a 
disordered one for ri > and r2 > 0; 2) an ordered phase in which 0i 7^ and 
02 = 0, with ri < and r2 > ^Jbo/aQri, 3) an ordered phase in which 0i = 
and 02 7^ 0, with r2 < and r2 < a/^o/^o^^i- The multicritical point corresponds 
to Ti = r2 = 0, and, according to standard terminology, is named bicritical in this 
case. The phase boundaries between the disordered phase and the ordered ones 
correspond to second-order transitions with symmetry breaking 0(4) — > 0(3), while 
the boundary r2 = a/^o/'^o'"!, Ti < 0, r2 < 0, correspond to a first-order transition. 
In the second case, there are four phases: 1) a disordered one for ri > and r2 > 0; 
2) an ordered phase in which 0i 7^ and 02 = 0, with ri < 0, bounded by a line 
depending on the quartic parameters; 3) an ordered phase in which 0i = and 
02 7^ 0, with r2 < 0, bounded by a line depending on the quartic parameters; 4) an 
ordered phase in which both 0i 7^ and 02 7^ 0, bounded by the two lines considered 
above. The multicritical point corresponds to ri = r2 = 0, and is called tetracritical. 
The phase boundaries always correspond to second-order transitions. 

The nature of the phase diagram depends on the values of the quartic param- 
eters. The multicritical point is bicritical for cq > 2^/aQb^ and cq + > 2y/aobo. 
In all other cases it is tetracritical. In the tetracritical case, it is of interest to 
understand the symmetry breaking observed at each transition. The transitions be- 
tween the disordered phase and the ordered ones always correspond to the breaking 
0(4) ^ 0(3), or more precisely to Z(2) ® 0(4) Z(2) ® 0(3). The additional 
Z(2) corresponds to reflections of the nonmagnetized field with respect to the plane 
orthogonal to the magnetized one: for instance, if 0i is magnetized in the first di- 
rection, then the symmetry (02)i — > — (02)i) (02)i ~^ (02)i {i > 2) is not broken. 
In the presence of fluctuations these transitions are expected to belong to the 0(4) 
universality class. The transition lines bounding region 4 may correspond to differ- 
ent symmetry-breaking patterns. If do < and —2y/aobo < cq + do < 2^/aQbo, in 
phase 4 0i and 02 are parallel. At the boundaries the Z(2)-symmetry we discussed 
above is broken and the symmetry gets reduced from Z(2) (g) 0(3) to 0(3). There- 
fore, these lines should correspond to Ising transitions. On the other hand, if rfo > 
and —2y/aobo < cq < 21/0060 in phase 4 0i and 02 are orthogonal. Therefore, at the 
transition one observes the symmetry-breaking pattern 0(3) — > 0(2) (more precisely 
Z(2) ® 0(3) Z(2) 0(2)). In this case, the critical behavior should belong to the 
0(3) universality class. 

Finally, let us discuss the phase diagram predicted by mean-field theory in the 
limit in which the effective breaking of U{1)a is small. If g parametrizes this breaking, 
then Co = 3uo+^vo+0{g), do = —3vo+0{g), and ^/a^ = ^{uo+vo/2)+0{g), so that 
co + do = 2\/aQbo + 0{g). Using the results reported above, we obtain that for vq > 0, 
i.e. do < 0, the phase diagram is bicritical or tetracritical depending on the sign of 
the corrections proportional to g. Instead, for vo < 0, i.e. do > 0, the phase diagram 
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is tetracritical. We recall that the case relevant for QCD should be that with fo > 0, 
because only in this case the symmetry breaking \U{2) ® U{2)]/U{1) U{2)/U{1) 
is realized for g = 

The results we have obtained above in the mean-field approximation may change 
when we take into account fluctuations. As discussed in Ref. |]10[, fluctuations change 
the phase boundaries that are now generic curves which meet tangentially at the 
multicritical point. It is also possible that some of the transitions become of first 
order. This should happen in the case of interest of QCD. If g is the anomalous 
breaking of the axial U{1) symmetry, we know that for g = Q the transition is of 
first order. Moreover, (7 = should be the multicritical point, even in the presence 
of fluctuations, since such a point has a larger symmetry group. Since a first-order 
transition is generally robust under perturbations, the transition should maintain its 
first-order nature also for sufficiently small values of \g\, therefore the transition lines 
ending at the multicritical point are expected to be of first order sufficiently close to 
it. Possible phase diagrams are sketched in Fig. 13. 
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Table 3: The coefficients 6^^-^, 6^^.^, and 6^-^, cf. Eqs. Q and Q. Those 

corresponding to values of i,j, k that are not reported are zero. 
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